Magnetotransport through graphene nanoribbons at high magnetic fields 
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We have investigated the magnetoresistance of lithographically prepared single-layer graphene 
nanoribbons in pulsed, perpendicular magnetic fields up to 60 T and performed corresponding 
transport simulations using a tight-binding model and several types of disorder. In experiment, 
at high carrier densities we observe Shubnikov-de Haas oscillations and the quantum Hall effect, 
while at low densities the oscillations disappear and an initially negative magnetoresistance becomes 
strongly positive at high magnetic fields. The strong resistance increase at very high fields and low 
carrier densities is tentatively ascribed to a field-induced insulating state in the bulk graphene leads. 
Comparing numerical results and experiment, we demonstrate that at least edge disorder and bulk 
short-range impurities are important in our samples. 

PACS numbers: 72.80.Vp, 73.43.Qt, 73.22.Pr 



I. INTRODUCTION 



For the application of graphene in nanoelectronics one 
has to understand the behavior of graphene nanostruc- 
tures, in particular graphene nanoribbons (GNRs). They 
were theoretically predicted to show either metalhc or 
insulating behavior around the charge neutrality point, 
depending on their crystallographic orientation. In ex- 
periment, however, GNRs always exhibit an insulating 
state close to the charge neutrality point (CNP^, which 
is dominated by disorder rather than a confinement- 
induced gap in the spectrum^ ^ . A clear proof of conduc- 
tance quantization only appeared very recently in ultra- 
clean suspended nanoribbon^. Furthermore, i n c lean 
zigzag edges, a magnetic state has been predictecP^, but 
so far it has remained elusive in transport experiments. 
At present, therefore, the behavior of GNRs is mainly 
governed by extrinsic defects rather than their intrinsic 
properties, and information on the nature of those defects 
is highly desired. 



In previous experiments, large disorder was attributed 
to cause strong localization effects which influence the 
magnetoconductance^. Poumirol et al. report a large 
positive magnetoconductance and explain this by simu- 
lations which take into account different types of disor- 
der. They affirm the qualitative behaviour, but the com- 
puted conductance remains larger than the experimental 
ones. Also, an unambiguous separation of bulk and edge 
disorder was not possibl^. Here, we present magneto- 
transport measurements on GNRs in magnetic fields of 
up to 60 T and corresponding tight-binding simulations 
with several types of realistic bulk and edge disorder. By 
considering the magnetoconductance close to the Dirac 
point and at high densities, we observe characteristic sig- 
natures of bulk and edge disorder and can disentangle 
their contributions to transport in GNRs. 




FIG. 1: (Color online) (a) Scanning electron microscope im- 
age of a typical sample. The length of the GNRs is 1 /xm, 
the width 70 nm. In the upper part of the image a palladium 
contact is visible, (b) Two-terminal resistance as a function 
of Vhg at T= 25 K and zero magnetic field, (c) Magnetoresis- 
tance trace at Vbg = —20 V, showing quantum Hall features 
at = 6, 10 and 14. 



II. EXPERIMENTAL DETAILS 

Single-layer graphene is deposited on a highly doped 
silicon wafer with a 300 nm thick Si02 layer by conven- 
tional exfoliation. The graphene nanoribbons were de- 
fined by electron-beam lithography and oxygen plasma 
reactive ion etching. For the transport measurements, 
palladium contacts were attached to the GNRs. A scan- 
ning electron micrograph of the sample discussed here is 
shown in Fig. [ija). The DC magnetotransport measure- 
ments with 10 mV DC bias were done in pulsed perpen- 
dicular magnetic fields at temperatures between 1.8 and 
125 K. Typical pulse durations were ranging from 100 to 
300 ms. During the pulse the current through the GNR 
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was converted to a voltage signal by a current-to-voltage 
amplifier and recorded by a high-speed oscilloscope and 
data recorder. In total two single-layer nanoribbons have 
been measured which show similar behavior. Here, we 
focus on data from one device. Figure [T]jb) shows the re- 
sistance R of the nanoribbon as a function of back-gate 
voltage Vbg at T= 25 K and zero magnetic field. The 
sharp peak at Vhg — Vcnp = —4.4 V indicates the charge 
neutrality point. After patterning, the hole mobility /x of 
the ribbons is about 590 cm /Vs at Vbg= -15 Figure 
[ijc) shows a magnetoresistance curve taken at high car- 
rier densitjffSi A quantum Hall plateau dX v = 6)^1 and 
Shubnikov-de Haas oscillations for = 10 and 14 are ob- 
served. Signatures of Hall states were already found in 
previous experimentfi^!. From the zero-field mobility and 
the condition /ii? ^ 1 we would not expect to observe 
quantum Hall features ai v = 14, at 13 T. This is already 
an indication that the high field changes the impact of 
disorder on transport in our sample. 



III. DENSITY AND TEMPERATURE 
DEPENDENCE 
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FIG. 2: (Color) (a) Magnetoresistance of the CNR for various 
temperatures at the charge neutrality point, (b) Magnetore- 
sistance for different gate voltages close to the CNP and (c) 
further away from the CNP at T= 25 K. The arrows and 
the numbers indicate the corresponding filling factors of 
the quantum Hall state, v= 2 and 6. (d) Conductance as a 
function of magnetic field for Vbg= -15.6 and -6.2 V. 



Let us now consider the density and temperature de- 
pendence of the magnetoresistance in more detail. First, 
we will focus on the transport properties at gate volt- 
ages close to the CNP. For all temperatures we tuned 
the backgate voltage such that the samples remained as 
close as possible to the CNP. In Fig. [2ja) , the magnetore- 
sistance is plotted for various temperatures ranging from 
1.8 to 125 K. For all temperatures a resistance decrease 
is observed for fields up to about 20 T, so that the rib- 
bon crosses over from a highly resistive state to a metallic 
regime. Subsequently, it is followed by a prominent resis- 
tance increase. The divergent form of the latter increase 
suggests that the nanoribbon approaches a field-induced 
insulating state. 

In order to better comprehend the observed behav- 
ior, we studied the magnetoresistance for different gate 
voltages ranging from -4.8 to -13.7 V at T= 25 K. As 
one can see in Fig. [2]jb), the observed divergence of 
the resistance at very high fields only appears for gate- 
voltages close to the CNP {\Vbg -Vcnp\ < 9 V). At 
higher densities [see Fig. [2]^c)], we observe weak localiza- 
tion at fields up to 1 T, a fairly constant resistance up 
to about 20 T, and then pronounced resistance oscilla- 
tions. These oscillations can be identified as Shubnikov- 
de Haas (SdH) oscillations, which can be assigned to 
Hall-plateau values of single-layer graphene {v = 2 and 
6). The capacitive coupling Cg of the nanoribbon to the 
back-gate, which strongly depends on the ribbon dimen- 
sions, was calculated using a finite-element model, yield- 
ing Cg = 576 aF/jim^ for a 70 nm wide CNR. Plotting 
the fan diagram of the minima of the SdH oscillations 
gives a coupling Cg of 560 aF//xm^, which matches the 
calculated value well. Therefore, the carrier density is 
estimated as n « 3.5 x 10^^ m^^ x [Vbg — Vcnp) and the 



Fermi-energy scales as Ep 69 meV x yJ\Vbg — Vcnp\i 
where Vbg and Vcnp are given in Volts. 

For easier comparison to the numerical calculations. 
Fig. [2]^d) shows the conductance G as a function of mag- 
netic field for two different carrier densities representative 
for the low- and high-carrier-density regime. The high- 
carrier-density conductance {Vbg= -15.6 V) shows the 
oscillating behavior as described before, the low-density 
trace {Vbg= -6.2 V) exhibits first a conductance increase 
followed by a conductance decrease. In the following, we 
discuss the observed behavior with the help of numerical 
simulations. 



IV. NUMERICAL TRANSPORT SIMULATIONS 

The experimental data in Fig. [2] will give us impor- 
tant insight into the nature of the defects relevant in 
our GNRs. Specifically, in this section we will focus on 
the visibility of the SdH oscillations, the positive mag- 
netoconductance at low carrier densities and fields up 
to about 20 T, and the rather high zero-field resistance 
at both low and high carrier densities. To this end, we 
have performed numerical magnetotransport simulations 
of (armchair) graphene nanoribbons with realistic sizes 
(L = 320 nm, W ~ 25 nm). Since Ohmic scaling is not 
applicable at those length scaleJ^^ we do not expect a 
full quantitative match between theory and experiment. 
However, the qualitative behavior will be well reproduced 
by the simulations since the system size is of the same or- 
der as the experimental samples. We used the well-known 
graphene tight-binding Hamiltonian in nearest neighbor 
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so-called electron- hole puddles, i. e., long range potential 
fluctuations due to charged impurities trapped beneath 
the graphene ribbon in the silicon-oxide substrate. Sec- 
ond, we also consider shorter-ranged impurity potentials, 
that can arise due to adsorbates, defects or charged im- 
purities. In both cases, we add Gaussian on-site poten- 
tials to the tight-binding Hamiltonian Q. For the pud- 
dles, we use Gaussians with a decay length of ~ 8.5 nm 
and a total height of ^ 80meV, which is comparable 
to the experimentally determined valueJ^. The impu- 
rities were modelled by Gaussians with a decay length 
of-0.44nnPl. 



FIG. 3: (Color online) Magnetoconductance of armchair 
GNRs (L = 320 nm, W ~ 25 nm) calculated numerically, us- 
ing tight-binding simulations^ and different disorder models, 
(a) Edge disorder (cf. text, inset: a close up of the ribbon 
edge with disorder), (b) Long-range Gaussian disorder (pud- 
dles, cf. text), (c) Short-range impurities. We used Gaussian 
disorder with a decay length of ~ 0.44 nm. The height of the 
individual Gaussian potentials is randomly distributed within 
the interval [—S, S] with 5 = 0.1 1 and the impurity density is 
p — 15%. (d) Edge disorder and short-range Gaussian disor- 
der. Here 5 = 0.09 1 and p = 8%. 

(n.n.) approximation, 

ij'n.n. 

where for finite magnetic field the corresponding hop- 
ping integral is given by tij = —texp[ie/h J^' dsA{x)], 
with constant t 2.7 eV and the vector potential A{x). 
The conductance was then computed using an adaptive 
recursive Green-function method, capable of treating ar- 
bitrarily shaped system^^^. 

To appropriately describe the experimental situation, 
we considered different types of disorder. Since the fab- 
rication process certainly leads to disordered edges, we 
took this into account also in the numerical simulations. 
To this end, we cut 'chunks' of about 4nm out of the 
graphene lattice at random positions close to the edge, 
which simulates the large-scale edge roughness that oc- 
curs due to e-beam resist roughness and the random na- 
ture of reactive ion etching. Additionally, we accounted 
for edge roughness on a smaller scale of a few lattice 
constants using a model introduced in Ref.'^ About 10 
percent of the edge atoms are randomly removed and sub- 
sequently dangling bonds are additionally removed. This 
procedure was repeated 5 times to yield an edge rough- 
ness of a few lattice constants. The numerical results, 
however, showed that both types of disorder yield simi- 
lar results. In the following, in the case of edge disorder, 
both mechanisms will always be included. 

In addition to the edge disorder, we studied two types 
of bulk potential disorder. On the one hand, we modeled 



In Fig. [Sj we present our numerical results for magne- 
totransport through disordered nanoribbons at relatively 
high {Ep w 226 meV) and lower {Ep w 92 meV) carrier 
densities, corresponding to the Fermi energies of the ex- 
perimental data in Fig. l2jd). First, we consider ribbons 
with edge disorder only [Fig. [sj^a)]. We find that while 
the zero-field conductance for low densities is compara- 
ble to the experiment, this is not the case for the high- 
density result. Upon increasing the field, the wavefunc- 
tions become more localized close to the edges. With- 
out bulk disorder, backscattering is strongly suppressed, 
so that calculations yield nearly perfect quantum Hall 
plateaus for all densities already at moderate fields, in 
contrast to the experimental findings. This means that 
edge disorder alone cannot explain the experiment. Con- 
sidering only long-range Gaussian disorder [panel (b)], 
we find that the puddles are rather effective scatterers at 
low density, while they affect G only little at high densi- 
ties. Simulations where only the short-range impurities 
are taken into account [panel (c)], show that indeed for 
strong enough scattering potentials, the zero-field con- 
ductance can be very close to the experimental data. 
However, such strong bulk disorder leads to backscat- 
tering even for very high magnetic field, so that at high 
carrier density no SdH oscillations can be observed. This 
implies that indeed a combination of bulk and edge dis- 
order is necessary to describe the high-field experiments. 
In panel (d), we show the results for ribbons with disor- 
dered edges and short-range bulk disorder. In this case, 
the experimental findings for low and moderate field are 
reproduced semi-quantitatively. For low density, we find 
a strong increase of G due to the formation of edge chan- 
nels, while clear SdH oscillations are obtained at higher 
densities. The zero-field conductance fits well with the 
experiment. In contrast, in simulations that additionally 
include the long-range puddles, the difference in the zero- 
field conductance for high and low densities is much too 
high, thus we conclude that puddles are not the domi- 
nant scatterers in our samples. We note that beyond our 
disorder model interaction effects may further influence 
the measured conductance. 
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V. HIGH FIELD INSULATING STATE AT LOW 
DENSITIES 

We now turn our attention to the sample proper- 
ties at high magnetic fields near the CNP. As shown 
in Fig. [2]^ a) , the resistance at low temperatures initially 
decreases with B and then diverges steeply by several 
orders of magnitude for i? > 20 T. While the initial neg- 
ative magnetoresistance at low densities is explained in 
the previous section by the formation of edge channels re- 
lated to the zero-energy Landau level (LL) in graphene, a 
crossover to a divergent resistance for i? > 20 T requires 
another transport mechanism. The zero-energy state in 
bulk graphene has been investigated by several research 
groups, and a strong increase in R at the CNP and in- 
tense magnetic fields has been observed, resulting in a 
B-dependent L L spli ttin ^^^ l -^^ l and eventually a strongly 
insulating state^SEB^ the exact nature of which is still 
under debat^^H 

Adopting a simple model involving the opening of a 
field-dependent spin gap—, we can fit the temperature 
dependence of R for T > 14 K in an Arrhenius plot for 
distinct magnetic-field values (inset of Fig. ^ . In Fig. [4j 
energy gaps, A, are extracted from linear fits to the Ar- 
rhenius plot. The gap A shows a linear dependence on B 
(Fig. [4]), consistent with spin splitting of the zero-energy 
LL, with the gyromagnetic factor g=1.73. However, an- 
other origin of the gap can also be considered. Following 
for example Ref.'S^, we can fit A cx C • (i? — Bc)°'^ with 
Be ~ 29 T and C « 11, see Fig. |4} suggesting a chiral 
symmetry breaking transition. Comparing these different 
models we conclude that both mechanisms are compati- 
ble with our data, but the exact nature of the gap cannot 
be determined experimentally. For lower temperatures 
(T < 7 K) , however, the resistance diverges strongly with 
B, and a simple activated behavior can no longer explain 
our data. This divergent behavior of R in our GNRs re- 
sembles a field-induced transition to a st rongly insulating 
state reported in bulk graphene at low TflSEl] jn cleaner 
samples the transition to the insulating state occured at 
significantly lower fields. 

Given the sample geometry displayed in Fig. TJa), we 
note that (bulk) graphene leads are attached to the GNR. 
Since our GNRs, after patterning, have lower mobility 
than the bulk graphene leads the field required for the 
B-induced insulating state is expected to be also higher. 
Therefore, the observed divergent R at very high B and 
low densities is tentatively attributed to the leads: when 
we apply high B-fields the leads become insulating and 
mask the electron transport in the GNR. 



VI. CONCLUSIONS 

In conclusion, we have performed transport experi- 
ments in graphene nanoribbons in pulsed high magnetic 
fields and corresponding transport simulations, based on 
a tight-binding model. This allows us to separate the 
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FIG. 4: (Color online) Energy gaps, A, extracted from the 
slope of the Arrhenius plot for T > 14 K (inset). The (red) 
dotted line fits the Zeeman splitting, A = {gjj,BB)/kB — 8.9 K, 
with the Bohr magneton /is, the Boltzmann constant fcs and 
a gyromagnetic factor of g — 1.73. The (blue) continuous line 
is a fit following Ref.'^, of. text. 



contributions of different disorder types to magnetotrans- 
port. At least a combination of edge disorder and short- 
range bulk impurities is needed to reproduce the experi- 
mental results semi-quantitatively. The short-range bulk 
disorder is responsible for the partial suppression of the 
quantum Hall effect, while the edge disorder, together 
with the bulk disorder, provides sufficient backscattering 
to explain the observed high resistance at zero field for 
all carrier densities. Additionally, we observe a magnetic- 
field-induced insulating state at very low densities, which 
presumably originates from the bulk graphene leads. 
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